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Abstract 

The a-function is a proposed quantity defined in even dimensions which has a 
monotonic behaviour along RG flows, related to the /3-functions via a gradient flow 
equation. We study the a-function for a general scalar theory in six dimensions, using 
the /3-functions up to three-loop order for both the MS and MOM schemes (the latter 
presented here for the first time at three loops). 
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1 Introduction 


There has been considerable recent interest over the possibility (hrst raised by Cardy 
[I]) of a four-dimensional generalisation of Zamolodchikov’s c-theorem |2] in two dimen¬ 
sions; namely a function a{g) of the couplings which has monotonic behaviour under 
renormalisation-group (RG) flow or which is dehned at hxed points such that auv —> 0. 
These two possibilities are referred to as the strong or weak a-theorem, respectively. A 
proof of the weak a-theorem has been proposed by Komargodski and Schwimmer |3] with 
further analysis in Ref. |1]. 

The strong a-theorem has been proved valid for small values of the couplings w, 
using Wess-Zumino consistency conditions for the response of the theory dehned on curved 
spacetime, and with a:-dependent couplings g\x), to a Weyl rescaling of the metric [7]. A 
function A is dehned which satishes the crucial equation 

diA = TjjP'^, ( 1 . 1 ) 

for a function Tjj which is dehned in terms of RG quantities and may in principle be 
computed perturbatively within the theory extended to curved spacetime and x-dependent 
g^ . Eq. fll.ip implies 

/i^A = P^^A = (1.2) 

a/i og^ 

where Gjj = T^ijy, thus verifying the strong a-theorem so long as Gjj is positive-dehnite, 
a property which holds at least for weak couplings in four dimensions. It is clear that if A 
satishes an equation of the form Eq. fll.ip then so does 

A' = A + gijf3^f3'^ (1.3) 

for any gjj (for a diherent T/j, of course). 

Further extensions of this general framework have been explored in Refs. [BHin]. We 
should mention explicitly here that for theories with a global symmetry, (3^ in these equa¬ 
tions should be replaced by a which is dehned, for instance, in Ref. [6]. However, it 
was shown in Refs. mm that the two quantities only begin to diher at three loops for 
theories in four dimensions. 

It was shown in Ref. [13] that equations of a similar form to the above may be derived 
(in a similar manner) for a renormalisable theory in six dimension^; though the dehnition 
of A and Tjj as renormalisation-group quantities is of course diherent. For instance, 

A = 6a + bi-^h + WiP^ (1.4) 

where Wj, like T/j, has a dehnition in terms of RG quantities, a is the /S-function corre¬ 
sponding to the six-dimensional Euler density 

fP — f: p ]D^ll^^2VlU2 p/i5/i6l^5^'6 (-\ K;'\ 

^ See Ref. [14] for an analogous extension to three dimensions; and Ref. [15] for a recent explicit 
construction of an a-function in three dimensions. See also Ref. m for attempts to derive a weak a- 
theorem in six and general d dimensions using the methods of Ref. [3] . 


1 







and bi and 63 are the /5-fnnctions corresponding respectively to 


ii = + 1^2 - IQ, U = -sisii', + jiiiA'2 - ±K2 + + iA'o, (1.6) 


SO"*-! I 4-^‘-2 -‘*■( 3 ) eooo 

where (following the notation of Ref. [T3] ) 

Ki = R\ K2 = RTs = 

?kA d _ TjnX 


K, = R^^R^^,^R^\ K, = R^^R^^^pR, 


AA/. 13 , K4 = R^^Rxf^R^"^, 

fiup 


( 1 . 7 ) 


We nse the notation A to avoid confnsion with the a of Ref. [13] which differs by a factor 
of 6 . In six dimensions, Gjj has recently been compnted to be negative dehnite at leading 
order for a mnltiflavor 0 ^ theory miiiH]. The six-dimensional case has also been considered 
in more general terms in Ref. [T9|. Onr pnrpose here is to extend the resnlts of Ref. [TTlfTS] 
to higher orders, again for the mnltiflavor 0^ theory. However, we shall do this by nsing the 
0-fnnctions together with Eq. fll.ll) to constrnct the qnantities A and T/j order by order 
(rather than by a direct pertnrbative compntation). We shall compnte the fnnction A np 
to 5 loop order in the standard MS renormalisation scheme, reqniring a knowledge of the 
three-loop MS /3-fnnction. We shall hnd that a solntion for A and Tjj is only possible if the 
0 -fnnction coefficients satisfy a set of consistency conditions, and we shall be able to show 
that these conditions are invariant nnder the conpling redehnitions which mnst relate any 
pair of renormalisation schemes. We illnstrate the redehnition process nsing the example 
of the MOM (momentnm snbtraction) scheme. We accordingly present for the hrst time 
the three-loop /5-fnnctions for MOM (the three-loop MS /5-fnnctions may be read off from 
the results presented in Ref. [20], although they were not written down explicitly ther^. 
We shall also provide full details of the three-loop calculation for the MS /9-functions, and 
then give a precise dehnition of the MOM scheme, explaining how the calculation may be 
adapted for this case. In the general case we shall be considering, the theory has a global 
symmetry; and just as in four dimensions, we shall hnd that at three loops the consistency 
conditions can only be satished if we replace {3^ by the quantity dehned in Ref. [ 6 ]. 

The layout of the paper is as follows: in Section 2 we present the one and two loop 
results for the /9-functions and also the lowest-order results for A and Tjj. In Section 3 we 
give an explanation of our computational methods and how they may be applied to the 
computation of the /9-functions in both the MS and MOM schemes, and then go on to list 
the results for the three-loop MS /9-function. In Section 4 we present our results for A^^'^ 
together with consistency conditions on the /9-function coefficients, which must be satished 
in any scheme in order for Eq. fll.ip to hold. In Section 5 we discuss the implementation 
of renormalisation scheme changes in general terms and then go on to focus on the case 
of the MOM scheme. We present our concluding remarks in Section 6 , and hnally the 
explicit three-loop MOM /9-function together with some calculational details are given in 
the Appendices. 

"‘Earlier three-loop results were presented in Ref. [HIES], but these do not determine the results for the 
general theory unambiguously, as we shall explain later. 
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2 One and two-loop results 

We consider the theory 

L = ( 2 . 1 ) 

involving a multiplet of helds 0® coupled via the tensor The one and two loop /3- 

functions are given by 

~ 9(la) + J29{1A)-, 

= C(2B)g(2B) + C{2C)9{2C) + C(2fe)5'(2b) + C{2c)9{2c) + C{2d)9{2d), (2.2) 

where the tensor structures are dehned by 

afa) = aft) = a‘”a'‘’”'a\ty 

9g, = 9””y"9;>‘". 9|2) = 9-’"V"9‘’"9""‘9’”’". (2.3) 

with also 

9;b)=9‘V“. 9gB,=9"<V Ujh) = 9™V’”''9ri%. (2.4) 

We also dehne 3-index quantities corresponding to the 2-index quantities of Eq. (12.dh by 

9^{i% = etc. (2.5) 

We have therefore given the same label to both a two-index and a three-index tensor, 
but we hope that it will always be clear from the context which is meant. For structures 
which are not three-fold symmetric, we list one symmetrisation but it is to be understood 
that it is accompanied by its two symmetrised partners. We shall also, wherever possible, 
suppress indices as we have done in Eq. fl2.2p . We display the tensor structures appearing 

the right, except for g( 2 d)-i which 


^ A 

9(2c) 9{2d) 

Table 1: Tensor structures appearing in the one- and two-loop /5-functions 
The coefficients in Eq. fl2.2p are given in MS by 

C(2S) — Ygj '^(2C) ~ ~ ~45 '^(2c) ~ 72’ *"(2ii) ~ ~ 2 ' 


in Eq. (12. 2 p in Table [U (in which the index i is always 
is completely symmetric in i, j and k). 



9{la) 9{1A) 



9{2B) 9{2C) 5 ( 26 ) 
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Here and elsewhere we suppress a factor of (ddvr^) ^ for each loop order. 

It was shown in Refs [T^IIH] that Eq. fll.ip was valid at leading order with 

yl® = - Igjf^,) = gtvglv, _ ig^Jjg””,), 

T?j = Gn = ( 2 . 7 ) 

and 

A = -5^. (2.8) 

As mentioned earlier, our dehnition of A differs from that of a in Ref. [13] by a factor 
of 6, introduced for convenience. As the notation implies, G\j would require a two-loop 
perturbative calculation using the methods of Ref. w, which was performed explicitly 
in Refs. imCH]; would correspondingly require a three-loop calculation, but its value 
was inferred by imposing Eq. (II.Ih . This is the technique we shall apply to obtain 

later in this paper. We should note that A has no explicit two-loop contributions; 
however there is a one-loop (free-held) contribution [T71ITB] El. At this point we have all the 
information required for our computation of A^^\ but we shall postpone this to Section 4 
where we shall explain the general method we shall use to compute both and A^^b 


3 Three-loop results 

In this section we shall give our explicit results for the three-loop /3-functions for the theo^ 
Eq. (12.ip . computed using MS. We start by describing the computational methods useqj, 
and we shall also take the opportunity to describe the MOM scheme (which will feature in 
later sections) and how to adapt our methods to obtain the MOM /3-functions. The reader 
uninterested in the technical details of the Feynman diagram calculations may skip to the 
paragraph containing Eq. (13.61) . 

We base our calculations on the original work of Refs. [2T1[22] , although our couplings 
are more general than those articles. However, one can always make contact with the results 
of Refs. [2I1I22] by setting g was the coupling constant and was a 

group valued object which is completely symmetric in its indices like For example, if 
0 * took values in the adjoint representation of SU{Nc) then (pA would be the corresponding 
totally symmetric rank 3 colour tensor of that group. From the point of view, however, 
of constructing the RG functions it is the evaluation of the underlying Feynman integrals 
which is required. For this aspect it is appropriate to focus for the moment on the basic 
((A Lagrangian 

L = \ + l4? . ( 3 . 1 ) 

®See Refs. [23l[2l] for earlier perturbative calculations on six-dimensional curved spacetime. 

® Similar methods have recently been applied to the computation of the four-loop MS /3-function for 
this theory pS] . 
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Previously this was considered in Refs. [2Tl|22l|26] to three loops and the wave function 
and coupling constant renormalisation constants were deduced from the divergences in the 
2- and 3-point functions. For the latter the pole structure was determined by exploiting 
specihc properties of the six dimensional spacetime. Briefly it was possible to nullify the 
external momentum of one of the legs of the 3-point functions to reduce the evaluation 
of the graph to a 2-point function. Such integrals are more straightforward to determine 
through knowledge of the 2-point evaluation. The point is that ordinarily the nullihcation 
of an external momentum leads to infrared problems. For instance, a propagator in a 
Feynman integral of the form leads immediately to infrared singularities in the 

associated four-dimensional integral where k is the loop momentum. In six dimensions, 
however, this is not the case. The corresponding propagator where such an infrared issue 
would appear in that dimension is Therefore, the nullihcation used in Ref. [22], 

which is a simple application of the infrared rearrangement technique, is perfectly valid in 
six dimensions for Eq. fl2.1l) and Eq. fl3.1l) . Moreover, this method is sufficient to determine 
the RG functions in the MS scheme. 

However, as our focus here will not be restricted to MS but will include the momentum 
subtraction (MOM) scheme we will have to determine the Feynman integrals contributing 
to the vertex function for a non-nullihed external momentum and to the hnite part in e 
(where d = 6 — 2 e). Ahead of the description of the computational tools we use it is appro¬ 
priate to recall the dehnition of the MOM scheme as this informs the integral evaluation. 
The particular MOM scheme we use is that developed in Ref. [27| for Quantum Chromo¬ 
dynamics. However, we note that it was used prior to that at three loops for Eq. fl2.1l) in 
Ref. [28| when the specihc coupling tensor corresponded to an SU{3) x SU{3) symmetry 
group. First, we recall that in the minimal subtraction scheme the renormalisation con¬ 
stants are dehned at a subtraction point in such a way that only the poles with respect to 
the regularizing parameter are included in the renormalisation constant. The MS scheme 
is a variant on this where a specihc hnite part is also absorbed into the renormalisation 
constants. This additional number, which is ln( 47 re“'’') where 7 is the Euler-Mascheroni 
constant, in ehect corresponds to a trivial rescaling of the coupling constant. The renormal¬ 
isation constants of the MOM scheme by contrast are dehned at a particular subtraction 
point in such a way that after renormalisation at that point there are no 0{g) contributions 
to the Green’s functions m In other words the Green’s functions are set to their tree 
values at the subtraction point [27]. The specihc subtraction point used in Ref. [27] for 
the 3-point vertex functions is that where the squared momenta of all the external legs 
are equal. Moreover, they are set equal to (—/i^) where /x is the mass scale introduced to 
ensure the coupling constant or tensor is dimensionless in general d dimensions. To achieve 
this for Eq. fl2.1l) requires evaluating all the 2-point function and vertex graphs to the hnite 
parts in e; the latter being a more involved computation than the former. While this is the 
canonical dehnition of the MOM scheme we note that one can have variations on it. For 
instance, it is acceptable to have a scheme where the renormalisation constants associated 
with 2-point functions are dehned in an MS way but the coupling-constant renormalisations 
are determined using the MOM dehnition or vice versa. While it is possible to study such 
hybrid schemes in order to consider diherent applications of our general formalism, we will 
concentrate here purely on the MS and MOM schemes. However, we will provide enough 
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details in the evaluation of the symmetric-point 3-point integrals to allow an interested 
reader to explore these hybrid schemes independently. In mentioning that a MOM scheme 
had been considered in Ref. [28] we need to clarify this in light of our discussion. The 
renormalisation of Ref. [2H] was at three loops in a MOM scheme but for a very specihc 
coupling tensor. The consequence of the choice of the SU{3) x SU{3) colour group is that 
there were no one or two loop vertex graphs to be evaluated at the symmetric point. At 
three loops due to the special symmetry properties there was only one graph to evaluate 
at the symmetric point, and since it was primitively divergent the evaluation was straight¬ 
forward [28|. Therefore, in considering the general cubic theory in six dimensions Eq. fl2.ip 
we are hlling in the gap in the computation of the lower loop integrals for the momentum 
subtraction scheme analysis. 

One of our aims is the study of the a-theorem up to third order in various RG schemes; 
it may appear a formidable task to actually calculate the MOM /3-function for Eq. fl2.ip 
at this order given the previous discussion. However, it is possible to determine it purely 
from the two loop renormalisation of the vertex functions in this scheme. To illustrate the 
process, we consider for the moment the simpler theory Eq. fl3.ip . postponing the general 
case of Eq. fl2.ll) to Sect. 5. After choosing a renormalisation scheme for a given theory, any 
RG quantity will be a function of the couplings in that scheme. However, the expressions 
in different schemes must be related and this is achieved by a conversion function. If 
we denote the coupling constant in one scheme by g and that in another by g then the 
conversion function Cg{g) is dehned by 


CM 


dg(g) 

Sg 


(3.2) 


where the bare coupling 

constant go 

is related to the two i 

by 




9o -- 

— p. Zg g P Zg g . 

g{g) has the form 


OO 


Kg) 

= S + 


(3.3) 

(3.4) 


n=l 


where the coefficients hn are related to the hnite parts of the Zg and Zg. Once the expansion 
has been established from the explicit renormalisation then the respective /3-functions are 
related by 

Kg) = [Kg)C9i.g)]g^g (S-S) 

where the mapping means that the coupling constant g is mapped back to g via the inverse 
of Eq. fl3.4p . Now returning to our problem of computing the MOM /3-function, if g is the 
MOM coupling constant and g is that in the MS scheme, it is clear from Eq. (13.51) that 
only hi and ^2 are required to hnd Kh) three loops . These may be derived from the 
hnite part of Zg up to two loops which in turn derives from the hnite parts of the one and 
two-loop vertex functions at the symmetric subtraction point. The coupling redehnition 
required to generate the three-loop MOM /3-function will be presented in Sect. 5, and the 
three-loop MOM /3-function itself is given in Appendix A. 
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We now turn to the algorithm we used to evaluate the two and three loop Feynman 
integrals to the requisite orders in e to determine the RG functions in MS and MOM. The 
method we use is to apply the Laporta algorithm [29] to the two and three loop graphs 
contributing to the 2-point and vertex functions. This method systematically integrates 
by parts all the graphs in such a way that they are algebraically reduced to a basic set 
of what is termed master integrals. Then the e-expansions of the latter are substituted 
to complete the computation. The masters have to be determined by direct methods or 
one which does not use integration by parts and this is the more demanding aspect of the 
calculation. The version of the Laporta algorithm which we used was Reduze [301131]. To 
handle the surrounding tedious algebra we used the symbolic manipulation language F ORM 
[32l[33]. The whole evaluation proceeded automatically by generating all the Feynman 
diagrams electronically with the Qgraf package [M]. While we have summarised what 
is now a standard procedure to carry out multiloop Feynman graph evaluation, the novel 
feature here is hnding a method to access the six dimensional master integrals. The most 
straightforward way to proceed is to realise that if the problem was in four dimensions 
then the corresponding masters are already known. For instance, the three loop 2-point 
functions were developed for the Mincer package [35] which was encoded in Form in 
Ref. [36]. Equally the two loop 3-point function masters were determined over a period 
of years in Refs. [37fHT] . While the algorithms which were developed to determine such 
four-dimensional masters could in principle be extended to six dimensions, in practice this 
would be tedious. Instead we have used the Tarasov method developed in Refs. [i2llT3] . 
This allows one to relate a scalar (master) integral in d-dimensions to integrals with the 
same propagator topology in (d -|- 2) dimensions. The caveat is that the integrals in the 
higher-dimensional case have increased powers of their propagators. However, the Laporta 
algorithm [29] can be applied to them in order to reduce them to the corresponding master 
in that dimension and integrals which involve masters where the number of propagators 
in the topology have been reduced. In summary a four dimensional master can be related 
to its unknown six-dimensional cousin plus already determined lower-level six-dimensional 
masters. Hence one can algebraically solve the system for the masters to the order in e 
which they are required for the 2-point and vertex functions. One feature which invariably 
arises in the use of integration by parts is the appearance of spurious poles. Consequently 
in the determination of the masters in six dimensions as well as the reduction of the 
Feynman graphs contributing to the RG functions, one sometimes has to evaluate the 
masters beyond the hnite part in e. In the appendix we record the expressions for the 
nontrivial one and two loop 3-point masters; though we stress that the expressions given 
there are for the pure integral. There has been no subtraction of subgraphs as was the 
case in the results presented for graphs in Refs. [201122] . This is primarily because in 
the automatic Form programmes to determine the RG functions we use the algorithm of 
Ref. [H] to implement the renormalisation automatically in our two schemes. Briefly this 
is achieved by performing the calculation in terms of the bare coupling constant or tensor. 
Then what would be termed counterterms are introduced automatically by rescaling the 
bare parameter to the renormalized one. The so-called constant of proportionality would 
ordinarily correspond to the coupling constant renormalisation constant which therefore 
appends the necessary counterterms. Finally, the dehnition of the renormalisation constant 
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is implemented at this last stage. At whatever loop order one is working to, the remaining 
undetermined counterterm is dehned according to whether the scheme is MS or MOM. 

We can now present our three-loop results for MS. At this order we require several new 
tensors. Those contributing to the one-particle irreducible (IPI) terms in the /^-function 
are 


ijk 

9{3e) 


Aim Alkm 
i/(26)i/22 5 


ijk 

9{3h) 


- 9{la)9^iA)9 9 


ijk 

9{3k) - 


„mil Alkm 

9{2c)922 ) 


ijk 

9{3n) 


— AlnAmq klm nq 

— y y y y{ 2 D)^ 


ijk 

9(3q) - 


„irs Jqps kpqr 
”(la)i/22 i/22 ) 


ijk inq jrps knrl nqls 

9{3t) = ^ T /22 £/22 £/22 , 


ijk 

9i3f) 

ijk 

9{3i) 

ijk 

9(31) 

ijk 

9(30) 

ijk 

9{3r) 

ijk 

9i3u) 


Imi jlkm ijk _ Apqjpr kqr 

i/(2fe)i/22 ) 9(3g) y i/(la)i/(la)’ 

Aim jlkm ijk _ Ain jmq klm nq 

i/(2c)i/22 ) i/(3m) — y y V ^(20)’ 


„Pq „rs Jpkr iqs ijk 

9(iA)9{iA)y22 y ■> y{ 3 p) 


y i/(lyl)i/22 i/22 ) 


Apq „jrs knps rnq ijk _ Um Jlkm 

y y y22 y(la)’ 9 ( 3 s) — 9{2d)922 ! 

kva jvsn ispl rlqn 

9 y22 922 922 , 


(3.6) 


where 


ijkl _ Ajm klm ij _ im mj /o 7 ', 

^22 —9 y 1 9{2D) ~ 9{lA)y(lA)i l'J-0 

and (/(la), y(iA) etc are dehned in Eqs. fl2.3p . fl2.4p . The tensors contributing to the anoma¬ 
lous dimension are 


y(3^) = 

9{3E) 

= y^F) 

— jpqjP’i 
~ 9 y(2c)’ 

— Amn mp pnj 

9 9(iA)9{la)-i 

ij 

9{3G') 

— n^P 

y(ia)y{iA)y ’ 

ij _ ipjqA‘1 

9(3H) ~ 922 y(2B)’ 

y(37) = 9"^A\ll). 

ij 

9(3J) 

— jrnn mp nq jpq ij _ ipjq pq 

— y y(iu)y(iyi)y ’ 9(3k) — y 22 y( 2 C)’ 


'Ij 

9{3L) 

injq nq 

— y22 y{2D)- 

(3.8) 


The notation for the two-index tensors in Eq. (13.8p matches the diagrams in Figure 7 
of Ref. [20], so that y^^^jg) corresponds to the tensor structure of Fig. 7(d), and so on. 

The notation for the three-index tensors in Eq. (13.6p similarly matches the 

diagrams in Figures 8 and 9 of Ref. [20]; and furthermore, the indices i,j, k are arranged to 
run anti-clockwise around the diagram, with i at the top. Finally, y( 3 ^)... y^^s^L), as dehned 
in Eq. fl2.5p in terms of ... ^(gx,), correspond to three-point diagrams with an insertion 
of the corresponding wave-function renormalisation diagram of Fig. 7 of Ref. [20] . 

As in Section 2, we write the three-loop /^-function as 

= C(3e)y(3e) + . . . + C(3u)g{3u) + C(3X))(/(3_D) -|- . . . + C(3L)y{3L)- (3.9) 

The coefficients in Eq. fl3.9p corresponding to anomalous dimension contributions are given 
by 


C(3D) -- 

„ — 121 
^{3H) 5184’ 


7 

864’ 


C(3E) 


71 

1728’ 


^{31) 96 24*5(3), 


C(3L) 


13 


C(3F) 

C(3J) 


10368’ '"(3G) C(3G') 

23 _ 103 

62208 ’ '^i3K) 7770, 


1 

108’ 


31104’ 


(3.10) 









where C(z) is the Riemann (^-function, and the remaining coefficients are given by 


C(3e) g) C(3j) 


^(3ff) ~ 16 ’ ~ 

47 

864’ 

C(3i) - 

47 

432’ 

'^(3i) - 

23 

288’ 

5 „ _ 

11 

C(3m) 

19 

C(3n) = 

11 

*^(3^) ~ 27’ *"(30 “ 

216’ 

324’ 

1728’ 

*"(30) ~ 1728 ’ *"(3p) ~ 

11 

144’ 

C(3g) = 

1 

16’ 

C(3r) = 

-i + C(3), 

C(3h = “if + iC(3), C(3t) = 

-1, 

C(3-!i) 

i-C(3) 


(3.11) 


We have compnted all the coefficients in Eqs. fl3.10|) . fl3.11|) explicitly and independently; 
and we have checked that we reproduce the wave function and /3-function results of Refs. 
[20] and |2T1I^ (in the latter case, after we specialize to the corresponding restriction on the 
group theory structure used there). Although in Ref. [20] the hnal /3-function results are 
given for two particular theories, the general results can be constructed from the individual 
diagrammatic results. This is largely the case for Ref. |211[22] too; however, the results for 
the pairs {c( 3 e), C( 3 /)}, {c( 3 fc), C( 3 p}, and {c( 3 d), c/s^)}, are presented together and cannot be 
separated. 


4 The a-function beyond leading order 

We now turn to the derivation of At two and three loops, we have avoided 

any diagrammatic computations using the methods of Ref. 13 E], as explored in the six¬ 
dimensional context in Refs. HU US]; instead we have proceeded to infer A by imposing 
Eq. fll.ip . Beyond leading order (corresponding to Eq. fl2.7p i we need to take into account 
potential higher order corrections to T/j. We suppress the details here but the calculation 
proceeds along similar lines to those presented in full in the four-dimensional case in Refs. 
[ISlIlG]; a similar method was used in the pioneering work of Wallace and Zia [17]. At 
next-to-leading order the general form of the a-function is given by 


= A 


1 + a^Af + a^Af + a^Af + a^Af + 


12 


(4.1) 
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Table 3: Next-to-leading-order metric terms 


where A is defined in Eq. fl2.8p . The individual contributions to depicted above in 
Table |2l are given by 

4(4) _ 4(4) _ ijk kij .( 4 ) _ ijk ijk 

—g 9{2d)^ ^2 —g fi'( 2 b)’ ^3 —g g{ 2 c)^ 

= s^^sgo). 4*’ = 9;4)9?4,9f;A). (4.2) 

with the tensor structures again defined in Eqs. fl2.3p . fl2.4p . fl2.5p . and (from Eq. fl2.2p l 

^il)ijk^{l)ijk ^ ^( 4 ) _ 1 ^( 4 ) ^ ^^( 4 ) ^ ^^( 4 ) _ 

Correspondingly the tensor Tjj in Eq. fll.ip is automatically symmetric at this order 
(so that Gij = Tjj) and may be written in the form 

4 

7 ^( 3 ) ^ '^ 43 ) /Tn( 3 )\ (AA'] 

ijk,lmn / \ a Jijk,lmn' 4 • J 

a=l 

where the individual structures which may arise are depicted in Table [3l Here the diagrams 
represent ta\Ta^)ijkjmn/3^^’^{dg)^"^'^ for a = 1.. .4. A cross denotes {dgY^^ and a diamond 
represents 

A careful analysis leads to a system of linear equations whose solution imposes a single 
consistency condition on the /5-function coefficients: 

6c( 2C) + C(2c) + C(2B) = 0. (4.5) 

This is satisfied by the MS coefficients as given by Eq. (12.6p . Similar integrability conditions 
(on three-loop /5-function coefficients) were found in Ref. [15]. In Eq. (14.ip . ai is arbitrary, 
reflecting the general freedom expressed in Eq. (II. 3p : and in particular for MS we have 

^2 8 ^1? ^3 ^ 2^1’ ^4 288 ^3 gg4 T ^^1* (^*^) 

Our methods therefore specify A up to the freedom expressed in Eq. fll.31) . The metric 
coefficients in Eq. fl4.4p are given by 

^2 ^ = —^A -|- ^oci, ^3 ^ T ^4 ^ = —|A -f- Oil, ^4 ^ = —6q;i, ( 4 - 7 ) 
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where once more A is defined in Eq. fl2.8p . The metric coefficients therefore also reflect 
the freedom expressed in Eq. fll.3l) : but there is an additional arbitrariness since only the 

/Q'j /Q\ 

combination ^ 3+^4 is determined. 

At the next order (corresponding to the three-loop /9-function) we need to face the 
possibility that {3^ in Eq. (inii should be replaced by a generalisation due to the 

invariance of the Lagrangian Eq. fl2.ip under 0{N) transformations of the real fields 0®. It 
was shown in Ref. [B] that in this situation, in the general case with couplings we have 

(3^ ^ = 13^ - {vgY (4.8) 

where v is an element of the Lie algebra of the symmetry group. In the case at hand, this 
corresponds to 

l^ijk Qijk ^ ^ijk _ ^l[igjk)l 

where v is an antisymmetric matrix. In principle v could be computed using similar 
methods to those described in Ref. | 6 ] and carried out explicitly in the four-dimensional 
case in Ref. |15]; but it is clear a priori that the relevant tensor structures are the same as 
those appearing in the three-loop anomalous dimension. Since most of those are symmetric, 
the only possible IPI contributions to v correspond to 

V = cYG)i9{3G) - g{3G')) (4.10) 

with g( 3 G), g{ 3 G') as defined in Eq. (13.81) . 

The solution of Eq. (II.ip leads to a complex system of linear equations whose solution 
imposes several consistency conditions on the 0 -function coefficients: 


C{3q) - C{3s) - 12C(3/) 
*^(3r) ‘3‘C(^3s) T 12C(3p) 

C( 3 e) ~ C(3g) — 24 C( 3 /j) — 144 C( 3 o) — 72Z 
C(3e) ~ C{ 3 g) ~ 6C(3j) -|- 6C(3fc) -|- 72Z 
^( 3 i) + 12c(3n) — C(3E) + 12C(3//) -|- 36 C( 3 i^) -|- 72 C( 3 L) 
“12(c(3p) -|- C(3o-)) — 72 (c( 3 t-) -|- C(3^)) -|- 6Z 
2C(3/i) -t- 6 c( 3 m) — 12C(3 ,i) -f 18 C( 3 o) + C( 3 D) + 12 C( 3 F) + 72 C( 3 J) 

+36c(3x) — 72c(3l) 

C(3e) — 1^(3/) -t- 6C(3fc) — 12C(3p 
^(3j) 4" 6c(3m) + 6C(3H) + 36c(3if) 
C{3h) — C(3i) + C(3;) — C(3£)) — 12C(3i?) -|- 12Z 


and 

12 C('3 c-) + 6(c(3 G) - C( 3 G')) + 12(c(3p) - C(3f^)) -h 72 (c( 3 r) - C( 3 ^)) 


—6C(2B), 

12C(2c), 

“3(C(2_B) + 2C(2c)), 
3c(2B) + 144c(2c) , 


12(2C(2c) + 2C(2c)C(2B) — C( 2 B)), 

^[1 + 24(c(2S) -C(2 c))], 

0 , 

12C(2S)C(2c), 

12C(2c)(C(2S) + 2C(2c)), 

(4.11) 

'^(3j)+ 6C(3m) -|- 12C(2c), (4-12) 
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where 


^ — C(3G) + C(3G') — C(3o) + ^C^SE) “ 20(3^) + 12C(3j). (4-13) 

We have included in our calculations potential contributions to the three-loop /9-functions 
dehned according to Eq. fl2.5l) in terms of one-particle-reducible (IPR) anomalous dimen¬ 
sion structures given by 

ij il Ij ij il Ij ij ik kj ij ik kj / a ^ a\ 

3{■ip) — 9{2B)9 {ia)^ 9{ia) — 9 {ia)9{2b)^ fi'{3r) “ 9{i^) — 9 {ia)9{2c)^ (4-14) 

using again tensors dehned earlier in Eq. fl2.4p . Such contributions cannot of course arise 
in MS but are potentially present in other schemes, in particular the MOM scheme which 
we shall be considering later as an example. (There are other potential IPR contributions 
in a general scheme, but we have included only the ones which are relevant for MOM.) In 
deriving Eqs. fid.lip . fl4.12p . we have imposed the two-loop consistency condition Eq. fl4.5p 
and also used the MS values of Eq. fl2.6p for all two-loop [3 function coefficients except C{ 2 b) 
and C( 2 c), since as we shall see later (in Eq. (I5.5p ). these are the only scheme-dependent 
values. The conditions in Eq. (14.lip are readily checked to be satished by the MS three- 
loop /9-function coefficients given in Eq. (I3.10p . (13.lip : however, Eq. (I4.12p may only be 
satished within MS by taking a non-zero value of namely 

(c0G))ms ^ “iMis- 

We shall assume this value from now on, though of course it would be reassuring to compute 
it directly along the lines of Ref. |T2] . 

Given that the consistency conditions are satished, we may solve Eq. fll.ip for the a- 
function. With the three-loop MS coefficients we hnd 


yl(5) 


where 






9 9{3p)i 


= 


~ 9 9{3r) - 


4f = 




4f = 

9'{[A)9{iA)9{iA)9'([A), 





= 

,(5) _ 




= 




4(51 _ 





Eof’d'’. (4.16) 


4(5) _ Jjk ijk _ ijk ijk _ ijk ijk 

4^9 —9 9{3M) ~ 9 9{3N) — 9 y(3n)’ 

4(5) ^ n^jk ijk ^ n^jkijk ^ Jjk ijk 

^12 9 5'(3e) 9 9 {‘if) 9 9{ig)i 

4(5) _ Jjkijk _ „ijk ijk _ „ijk ijk _ ijk ijk 

^13 ~ 9 9{3D) ~ 9 9{3h) — 9 9{3k) — 9 9{ 31)1 

4(5) _ ri^jk ijk _ Jjk ijk 

^14 9 9{3E) 9 9{3i)i 

4(5) _ „ijk ijk _ f.ijk ijk 

^15 9 9{3o) 9 9{3G)-i 

4(5) _ ijk ijk 
^16 ~ 9 9{3F)^ 

12 


( 4 . 17 ) 
























with the tensor structures etc as defined in Eq. fl3.6l) . The invariants are 

depicted in Table HI The coefficients in Eq. fl4.16p are given by 


(5) _ 
'1 ~ 

- 

V64 

sC(3)) A - 

4*^1 ’ 

Uq 

(5) 

'2 ~ 

( 29 

V 48 

+ iC(3)) A 

+ Cll, 

Uio 

(5) 

'3 ~ 

(1- 

iC(3)) A, 


“ii 

(5) _ 
'4 “ 

A + ^ do 

82944 ' ^ 144^2 

+ Ti4®3’ 

“l2 

(5) _ 

'5 ~ 

5 \ 11 

41472^ 864^! 

+ ^d2 + ^d3, 

“l3 

(5) _ 
'6 ~ 

29 \ 
2304^ 

- ^ai + 

M®3’ 

Ui4 

(5) _ 
7 ~ 

72®2’ 



“l5 

(5) _ 
'8 ~ 

-|A, 



“l6 


1152'^ 36^1 144®! gCl2 g«3, 

OQ -.I ~ 

~ 3®! ®3> 

” 1 ^"^ ^®1 ~ 3 ^ 3 ) 

+ !«! + di, 

” 32 "^ ~ 6 ®! ~ 3 ® 1 ; 

~^^A — |cri — ^di + d2, 

3456'^ + 72®1 ^®1 ~ 3^25 

~ 2304 144®! + 72®!’ (4-18) 


where A is again defined in Eq. (I2.8p and we explicitly display the freedom as expressed 
in Eq. (II.3p 

using the results of Eqs. fl2.2p . fl2.6p . The terms with di-d 3 correspond to the terms with 
in Table [3] (note that T^\ have the same effect at this order, as may be 
observed in Eq. fl4.7p b 
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Table 5: Next-to-next-to-leading-order metric terms 


At this loop order the tensor T/j is not inevitably symmetric. Analogously to the 
previous order, the metric may be expressed in the form 

30 

= ( 4 - 20 ) 

a=l 

The structures which may arise are depicted in Table El using a similar convention to 
the previous order. As before, a cross denotes {dgY^^ and a diamond represents 
We see immediately that ... t[^ are individually symmetric; the remaining diagrams 
are grouped in pairs whose coefficients should be equal for symmetry. Solving Eq. (Ea 
certainly does not guarantee that etc. However it turns out that we 

can impose symmetry on provided the additional condition 

C(3m) - 2C(3n) + |c(3£;) - C(3H) “ 12C(3L) - Z = ^{C(2B) “ 2C(2c)) (4-21) 
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is satisfied. The values in Eq. flS.lOp . fl3.1ip do indeed satisfy this condition and we obtain 
a symmetric metric with the values 


4") = iA-2ai, 

4 ^ ~ 24(52, 

4 ^ ~ —24'^ llCTl + 48(52, 

4 ^ ‘^*^1 24(52, 

4"^^ = + |(ni - 2(52, 

4*^^ = - 4d2, 

4"^ = -iA-2(52, 


^(4) _ _ n I 1 ~ 

^11 1728^ 72 3*^2, 

^(4) ^ J4_x I 4 ~ _ .x(4) 

^12 432^ ' 3^2 ^13 ; 

^(4) _ J_\ I 

^14 192^ ' 30 ^ 2 , 


^15 

_ 299 

1728 

A tiQ 

^''17 

+ 

+ §0^2, 



y-{4) 

y-(4) 







(-17 

~ ''18 ; 







'-19 

— _ 
~ '-20 ~ 

59 X J- 

72^4 + 

2.(3) 

3I4 

- 12((‘ie’ 

+ 47^ 

+ Tl"i 

+ 6Q!2 

^21 

— TT _ 
~ '-22 ~ 

115 \ 

288^ ' 

gtti — 

8(52, 




'-23 

— _ 
~ '-24 ~ 

73 \ 2 

48^ 3 

+ 

12(4? + 4?) - 

- i«i - 

12(52, 

^25 

— TT _ 
~ '-26 ~ 

101 \ 
288^ 

12®1 ~ 

■ 6(^2, 




'-27 

- - 
~ ''28 ~ 

373 \ 
1728"^ 

1U(3) 

36 4 

+ 2(4*6^ 

+ (‘l7 ) 

— f2^1 

— d2. 

^29 

— _ 
~ ''30 ~ 

+ 

1U(3) 

36 4 

- 2(4? 

+ 4?) 


2(52. 


The values of 4^^ 4?^ remain arbitrary, in a similar fashion to the previous 

order where in Eq. (14.7p only ^ 3+^4 was determined. Before imposing symmetry, the 
Tjf coefficients would also display the freedom Eq. fll.3p as expressed at this order in 
Eq. fl4.19p : however, the general redehnition Eq. fl4.19p is not compatible with symmetry 
of Tjf. Nevertheless, as may be seen in Eq. fl4.22p . there is still a residual two-parameter 
freedom expressed in (Ui, (52, corresponding to choosing 

(5i = —12(52 — d3 = 2(52 T (4.23) 

It is worth remarking that the freedom in ai is in general only preserved in the symmetric 
case providing Eq. fl4.5p is satished and is therefore somewhat non-trivial. We note that 
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in the four-dimensional case, the requirement of symmetry of Tjj was more restrictive and 
was not possible within MS |15]. 


5 Scheme changes 

In this section we shall turn to a fuller discussion of scheme changes such as that from MS 
to MOM. As we mentioned in Sect. 3, we have obtained the three-loop MOM /3-function 
by implementing the appropriate scheme change, avoiding a separate three-loop Feynman 
diagram calculation for MOM. Here we wish to consider the effect of more general scheme 
changes, in order to demonstrate the scheme-invariance of the consistency conditions on 
the /3-function coefficients, Eqs. fld.lip : we shall therefore give our results in full generality. 
We now rewrite the coupling redehnition of Eq. fl3.4p . which implements the change of 
scheme, in the form 

( 5 . 1 ) 

returning to the general couplings of Eq. fl2.ip and a general scheme change. The effects 
of Eq. fIS.ip may be computed from the generalisation of Eq. fl3.5p . 

(where /3 • ^ = /3^^™ ) which to lowest order may be written 

5/3*A = ^ ^ (5.3) 

dg dg 

The effect of a one-loop change 

Sg = Sig{ia) + S2g{iA) (5.4) 

on the two-loop /3-functions is easily computed as 

<5c(2s) = <^C(2c) = |A, A = 5i -|-1262 (5-5) 

where C( 2 S) and C( 2 c) are dehned in Eq. fl2.2p . It is readily checked that the consistency 
condition Eq. fl4.5l) is invariant under Eq. fl5.5p . as expected. 

At three loops we consider redehnitions 

^g = ^ig{2b) + e25'(2c) + ^3g{2d) + Q5'(2e) + ^3g{2f) + ^&g{2B) + ^7g(2C) + e85'(2D), (5.6) 

where 

9(2,) = g-'-gTiSaS). gg/ = gg/giL). ggf, = 9‘®9 (m)9S) (5.7) 
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Table 6: Two-loop IPR structures arising from coupling redefinitions 


as depicted in Table |6l and the remaining tensor structures are defined in Eqs. (12.3^ : and 
we also need to consider the effect at this order of the lower-order redehnitions given by 
Eqs. (El, (E5D . 

In general, in addition to modifying the coefficients already present in the MS /5-function 
as defined in Eq. (I3.9p (which correspond to IPI diagrams, or IPI wave-function renormali- 
sation diagrams attached as in Eq. fl2.5p l these redefinitions will generate tensor structures 
corresponding to IPR diagrams given by 

ijk il klj ijk jl Uk ijk U jlk ijk U Ijk 

9(!ia) = 9ilA)9^2iy ^?(3/3) = 9liA)9{2cy 9 ' 3 -y) = ^?(lA)^?(2c) > = 9(23)9 

ijk _ il Ijk ijk _ il jm Imk ijk _ il Ijk ijk _ il jm Imk 

9(3e) — 9{2C)9(ia)^ 9(30 — 9{lA)9(iA)9{la)y 9(37j) — 9{2D)9(ia)^ 9 ( 3 A) — 9 (IA)9 (23)9 , 


Jjk _ Am jl Imk ijk _ il jm Imk ijk _ il Am mjk 

9(3X) 9(iA)9(2c)9 ) 9(3^) 9(ia)9(2d) 9 i 9(3^) 9 (ia)9(2D)9 , 


(6.8) 


the tensor structures on the right-hand sides again defined in Eqs. fl2.3p . These IPR 
structures are depicted in Table. [71 we denote the coefficients of these tensor structures in 
the /5-function by csq, ... cs^, respectively. 

We obtain using Eq. fl5.2p (expanding now to 2nd order in (5i, 82 where necessary) 
^*^( 39 ) ~ ~2ei — 2 c( 2 b) 8 i + 81 , 

8C(3h) = |ei — £2 ~ C(2c)^l ~ ‘^C( 2 b )82 — 1 ^ 1 , 

8C(3i) = |ei — 262 + 265 — 2C(2c)5i — 4C(2fe)52 ~ ~ ‘^ 8382 , 

~ ~2e2 -|- 260 — 2 c(2c)5i + 2c(23)8i — — 451^2, 

8c(3k) = |ei -|- 262 + 2C(2 c)(5i — 4 c(26)(52, 

^C(30 = + £2 + C(2c)^l ~ 2c(2f,)52, 

8 c( 3 m) = ■|£2 + 2£7 — 4C(2c)^2 + 2C(2C)^1 ~ '^{\ 8 l 82 + 4(^2), 

8C(3n) = ^£2 + 268 ~ 4C(2c)52 ~ \ 8\82 — 78^, 

^C(3o) = ^£2 + £4 ~ 4 :C( 2 c )82 — \8182 — 5^2, 

<^C(3p) = |£3 — ‘iC( 2 d) 82 , ^C(3q) = “£3 ” C(2d)5i, 

8 C( 3 r) = —263 — 2 C( 2 d) 5 l, ^£( 3 ^) = 63 -|- C( 2 d) 8 l, (5-9) 


and 
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Table 7: Three-loop IPR structures arising from coupling redefinitions 


^C(3£)) = —— 2e6 -|- 40(26)^2 ~ 2 c(2b)5i -|- 4(5i52 + 
<5c(3e) = “1^1 ~ 266 + 2c(2b)52 ~ 2 c(2B)5i -I- -|- 4(5i(52, 

<^C(3i?) = —1^2 + 1^6 + 2C(2 c)^2 ~ 2C(2S)<^2 + 4^2 + 

<5c(3G) = i5c(3G') = —1^2 + 1^6 — ^7 ~ 1^5 

+ 2C(2c) 52 — 40(25)^2 — C(2C)^1 + |<5i52 -|- 86 ^, 
^C(3H) = “1^6 ~ 267 + 4C(2S)<^2 “ 2C(2C)<^1, 

<^C(3/) = —|e3 -f 2C(2d)52, 

^C(3J) = —|e4 + 1^7 — 2C(2C)^2 “ |<52) 

^C(3L) = ^67 — ^eg — 4 c(2C)<^2 “ 1^25 

and also for the IPR structures 


(5.10) 
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2 e 5 + 26162, ^C[3i3) = | e 5 — \6162, 6c(3^j = 2 e 4 + | e 5 — | 5 i 52 ~ 262, 

■2e5 + 26162, 5 C( 3 e) = ^€5 — ^6162, ^C(3ij) = —€4 + 62, 

1^5 ~ |<5i52, 5c(3k) = —4e4 + 4(52, '5 c(3a) = |e4 — |(52, 

^C(3^) = ^€4 — ^ 62 , 6c[3y) = leg — 1^2) 

<5c(3p) = —^€5 + ^66 — 268 ~ C{ 2 B )^2 + 4(52 + ^'^1*^2, 

^C(3o-) = ~ ^^6 ~ 2e8 + C( 2 B )62 + 2^2 + 

<^C(3r) = —1^4 + ~ C( 2 C )^2 “ |^2 

'^'^(3x) ~ “1^4 — + 1^8 + C{2C)^2 — |<^ 2 - (^-H) 

One may now verify that the consistency conditions Eqs. fl4.1ip are invariant under the 
coupling redehnitions of Eq. (15.dh : and so is Eq. (I4.12p . provided we assume that is 
scheme-independent, as is natural for a quantity appearing for the hrst time at this loop 
order. Finally, it is interesting to note that the condition required for symmetry of t)j\ 
Eq. fl4.2ip . is also invariant under these transformations so that a symmetric T/j can be 
obtained in any renormalisation scheme at this order. 

The coefficients C( 3 Q,)-C( 3 y) corresponding to IPR diagrams, and C( 3 p)-C( 3 ^) correspond¬ 
ing to IPR contributions to the anomalous dimension, are of course zero in MS. It would 
be natural to restrict ourselves to transformations which take us to other schemes with the 
same property; which indeed one might expect to be shared by any well-dehned diagram¬ 
matic renormalisation scheme. Now it is clear from Eq. (15.lip that any scheme change 
such that 

64 = 6 \, 65 = ^ 1 ^ 2 , ^8 = 1^2) (5.12) 

will preserve the vanishing of C{ 3 a)^C(^v), but will then inevitably give 

<^C(3p) = — 6 c { 3 ^) = g^ee + 62 — C(2b)<^2, 

^*^(3x) T 2 ^^ 6^2 C(2(7)^2- (5.13) 


^C(3q,) 

<5C(36) = - 


The values of C( 3 p), C( 3 o-) are opposite in sign, but not manifestly zero; likewise C( 3 r), C( 3 ^). 
Of course one could decide to choose eg, ey so as to ensure the vanishing of C( 3 p)-C( 3 ^) in 
Eq. fl5.13ll . but this seems rather restrictive and would not necessarily correspond to a 
natural diagrammatic renormalisation prescription. We shall therefore approach the issue 
from a different angle and consider in some detail the example of MOM, which certainly 
is dehned by a diagrammatic prescription as described in detail in Sect 3. 

The MOM and MS results hrst part company at two loops. The two-loop MOM /5- 
function coefficients may readily be computed directly using the methods described in 
Sect. 3. They are given by taking in Eq. fl2.2p 


MOM 

C( 2 B) 


1 _2_ 2 I 

36 81^ 27 




MOM 

C(2c) 


+ 


(5.14) 
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(5.15) 


where 'ip{z) is the Euler ■^-function dehned by 

the other two-loop /3-function coefficients remaining unchanged. Alternatively, we may 
compute the two-loop MOM /3-function coefficients by effecting the appropriate scheme 
change as described above. Comparing Eqs. fl2.6p . (15.141) . we simply require to take in 
Eq. O 

^ = 1 + Wil) (5-16) 

to effect the change to the MOM values. Of course this does not specify hi, 62 in Eq. fl5.5p 
uniquely. However, calculating the coupling redehnition required for the change from MS 
to MOM as described in Section 3, we hnd that the change to MOM corresponds to taking 

hi = | + Avr2-|i/i'(l), h 2 = -|, (5.17) 

and it is easy to check that Eqs. (15.5p . (I5.17P are compatible with the difference between 
Eqs. (12. 6 p and (I5.14p . 

The transformations required to change scheme from MS to MOM at the three-loop 
level may similarly be calculated using the methods described in Sect. 3. They are given 
by taking in Eq. (15. 6 p 

O = i + (1) , 

0 = 1- |C(3) - ^ ln(3)^\/37r 

- A ln( 3 )\/ 37 r - + 3^2 (f) - 6 S 2 (f) 

- 5s 3 (f) -f 4 s 3 (f) -F ^-0' (i) , 

O = -fi, O = tH, (5.18) 

together with 64 , 65 and eg as given in Eqs. (I5.12p . with hi 2 as dehned in Eq. (15.171) . 
The explicit MOM results for the three-loop /3-function coefficients obtained by combining 
Eqs. (I3.10p . (13.lip . (15.9p . (I5.10p are somewhat lengthy and are postponed to Appendix A. 

Unfortunately the values of eg, 67 shown in Eq. (I5.18p . do not correspond to the van¬ 
ishing of hc( 3 p), hc( 3 o-), hc( 3 ,-), hc( 3 ^) in Eq. (I5.13p . The scheme transformation therefore 
predicts non-vanishing MOM /3-function contributions from these IPR anomalous dimen¬ 
sion structures, which seems somewhat counter-intuitive. Indeed, after a careful direct 
calculation using the standard dehnition of MOM given in Sect 3, and taking account of 
the fact that the relation between /3-function coefficients and renormalisation constants is 
less trivial in MOM than in MS, we obtain = 0. It seems likely that the 

same applies to We therefore have an apparent inconsistency between the 

MOM values of C( 3 p)... C( 3 ^) obtained by the coupling redefinition process from MS, and 
those obtained by direct calculation. We have checked the anomalous dimension coeffi¬ 
cients obtained by coupling redefinition as given in Eq. (lA.lD by a 


£2 = 


703 


1728 


41 2 I ^ 

972 ^ 648 
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direct three-loop diagrammatic computation; it therefore appears likely that the discrep¬ 
ancy only affects the particular contriubtions C( 3 p)... C( 3 ^) corresponding to IPR anomalous 
dimension contributions. We also see from Eq. flb.lOh that this check of con- 

hrms the value of 67 and therefore hxes Sc( 3 t-) = —^C( 3 ^) 7 ^ 0 in Eq. (I5.13p . This removes 
the possibility that there might be a different choice of ei... 67 in Eq. flS.lSp which would 
correctly reproduce all the directly-computed MOM coefficients including vanishing values 
for C( 3 p)... C( 3 ^). It seems that one potential resolution of this problem lies in the use of 
a hybrid MOM scheme as alluded to briefly in Sect. 3, in which the wave-function renor¬ 
malisation constant is adjusted to give MOM values of C( 3 p)... C( 3 r) in agreement with the 
coupling redefinitions (without altering the values of any other /d-function coefficients). In 
six dimensions the issue may appear to be simply a technicality, but it seems probable that 
similar features arise in the four-dimensional case, which is of more practical interest. We 
propose to return to the subject in a subsequent article where we shall give full details of 
the MOM calculations reported here and show how a hybrid scheme can resolve the ap¬ 
parent inconsistencies. We shall also show how our results extend to the four dimensional 
case. Furthermore, it seems conceivable that a similar adjustment of the wave function 
renormalisation constant may be required in other schemes in order to match the results 
obtained by coupling redefinition to those obtained by direct computation, at least for IPR 
anomalous dimension contributions. 

We may now finally compute the MOM values of the coefficients in In order 

for this to be possible we know that the MOM coefficients should satisfy the appropriate 
consistency conditions derived in Sect. 4. We have already remarked that the two-loop 
MOM coefficients satisfy Eq. fl4.5p . Consequently, starting with A^^\ the MOM values of 
02 ... 05 in Eq. fl4.ip may be derived by again solving Eq. fll.ip . but with the MOM values 
as in Eq. fl5.14p : or more easily using the fact that A transforms as a scalar under coupling 
redefinitions, 

A{g) = A{g). (5.19) 

Using Eqs. fl5.4p . fl5.5p . and fl2.7p . we find for a general one-loop redefinition a leading-order 
change in A given by 


6A^^^ = A[A(kl2 - \A3) - 

and consequently for MOM using the values in Eq. fl5.17p we may take 


(5.20) 


^MOM 

“2 


= a2 + A = ^ + A7r2-|^'(i), 


a 


MOM 


= 03 - iA = A _ ^772 ^ (5.21) 


27 


with 04 and 05 unchanged. Clearly other choices are possible with corresponding adjust¬ 
ments of the value of ai in Eq. (14.ip . 

At the next order, we have mentioned already that the three-loop MOM coefficients in 
Eqs. (lA.ip . flA.2p satisfy all the consistency conditions in Eqs. (14.lip and (I4.12p . provided 
we take the non-zero MOM values of C( 3 p)... C( 3 ^) implied by Eq. (I5.13P (and listed in 
Appendix A), and assume that is scheme-independent, as is natural for a quantity 
appearing for the first time at this loop order. We have suggested above that these non-zero 
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values correspond to a hybrid MOM schem^. The MOM values of «(P • • • “pe) in Eq. (14T^ 
may then most easily be derived using Eq. flS.lQp and the explicit transformations given 
by taking Eqs. fl5.17p . flS.lSp in Eqs. fl5.4p . fl5.6p respectively. Again, one may also solve 
the equations using the MOM values of the /3-function coefficients as given in Eqs. fl5.14p . 
flA.ip . flA.2p : this will yield the same results, up to the freedom expressed in Eq. fl4.19p . 


6 Conclusions 

We have shown that, as in four dimensions, the gradient-flow equation Eq. fll.ip imposes 
constraints on the /3-function coefficients, and we have shown that these constraints are 
satisfied by the explicit results as computed for the MS and MOM schemes up to three- 
loop order. We have demonstrated that the tensor Tjj which appears in Eq. fll.ip may 
be chosen as symmetric up to this order. We have also shown that for a general scalar 
theory with an 0{N) global invariance, the /3-functions on the right-hand side of Eq. fll.ip 
must be replaced at three-loop order in the MS scheme by the generalised “S”-functions, 
as has also been observed in four dimensions. It would be useful to extend the analysis of 
Ref. [19] along the lines of Ref. |15] in order to understand the issues of theories with a 
global invariance further. This would have the beneht of enabling an explicit calculation 
of the “n” term in B (Eq. fl4.15l) ) and would also allow an understanding of its scheme 
dependence. 

Finally our analysis of scheme dependence has raised issues concerning the relation of 
MS and MOM; specihcally, the MOM values obtained for certain /3-function coefficients 
corresponding to IPR contributions to the anomalous dimension are different depending 
on whether they are obtained by direct calculation within the standard MOM scheme, or 
by coupling redehnition from MS. We shall discuss this issue further in a subsequent article 
where we shall show that the apparent discrepancy can be avoided by using a hybrid MOM 
scheme; we shall also explore similar issues in four-dimensional theories. 
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A Three-loop MOM results 

For the MOM scheme the three-loop /3-function coefficients were computed by constructing 
the appropriate coupling redehnitions as described in Sect. 3. Using the same notation as 

Ht is however interesting to observe that if we take C( 3 p) ... C( 3 ^) to be zero, as obtained by direct 
calculation in the standard MOM scheme, the MOM coefficients satisfy Eq. (14.1211 if we take = 0. 
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for the MS scheme, they are given by 
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MOM 
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MOM 

C(3G) 
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^MOM 
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C(3J) 


Cf 


432 486 ^ 2187 ^ 324^^ V3/ 729^" ^3v' ^ 243^" ^3/ ’ 
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M3G') 432 ^ 5832^ 3888^ Uh 
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- 1 ^ 2 ( 1 ) + S2(f) + § 53 ( 1 ) - is3(f) - 


2 
103 

31104’ 


MOM 

C(3L) 


85 

15552 ’ 


(A.l) 


for the anomalous dimension contributions, and 
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MOM 
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MOM 

^(3i) 
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^MOM 
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^MOM 
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C(3p) 
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162 108^ ''37’ 
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-A_L7r2 + 

288 324 ^ 216^ ''37’ 
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324 '' ^ 216 *^ '' 37 ’ 


288 


= i - K(3) - ^ 77 ^ + ^ ln( 3 ) 2\/377 - A ln( 3)\/377 - + S 2 (f) 


^MOM 

C(3g) 


- 2 s 2 (f) - 3^3(1) + 353(1) + 

= “i + iC(3) + - TIT hi( 3 ) 2\/377 + A ln( 3)\/377 + ^\/ 377 ^ - 3 s 2 (f) 


+ 6s2(f) + 5s3(f)-4s3(f)-2i^'(i). 


MOM 

C(3r) 


: -f| + 2C(3) + §77^ - ^ ln( 3)^\/377 + i ln( 3)\/377 + y|| 4 \/ 377 ^ - 6 s 2 (f) 
+ 12S2(f) + 1053(f) - 8 S 3 (f)- 7 ^'(i), 
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MOM 

(3s) 


t77^ + TIT ln(3)^\/377 - A ln(3)\/377 - oUgV^TT^ + 3 S 2 
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553 (f) + 453 ( 1 ) + 


(A. 2 ) 


for the IPI contributions. In Eqs. flA.ip . flA.2p . we define 
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Table 8 : One- and two-loop master integrals 


where ljin{z) is the polylogarithm function. The values of the coefficients C( 3 e), 

C( 3 /), C( 3 t), and C{^u) are identical in the two schemes. As mentioned in the main text, 
we have independently computed the values of initio and verihed 

that we obtain the same values. Also, the coefficients corresponding to the IPR anomalous 
dimension contributions obtained by coupling redehnition are 

C(3t) = ~C{3x) 


23 


10368 ’ 
61 


41472 • 


(A.4) 


B Master integrals 

In this appendix we record the explicit values of the various one and two loop six-dimensional 
master integrals which were needed to perform the renormalisation in the MOM scheme. 
Their derivation is based on the values of the corresponding master integrals in four di¬ 
mensions, which were given in Ref. [H], and are reproduced above in Tabled where Ad 43 ^ 
involves the square of a propagator, denoted by a “dot”. Using the same notation as 
Ref. IH] to denote the various graphs, the one-loop triangle master integral is 
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( 1 ) 

31 


13432 
2e 2 27 9 


23 


7r^ + 4s3(^) - 



35773 

324\/3 


ln( 3 )^ 7 r 

12\/3 


+ 0 ( 6 ^) . 


(B.l) 
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At two loops we have 
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^ ^ ^ 2^ 


-120^'(i) + 627r2 + 3915] 
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+ 7776s2(f) - 15552s2(f) + 10368s3(f) + 9667r2 - 4320C(3) 
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61 


15552 

1 25 205 7 2 1 ,//l^ n' 

12e2 ^ 48e ^ 96 ^ 648^ 27^ ^ 




1 

4^ 


+ 


'59 28 

^ - C( 3 ) - + 6.2(f) - 12.2(f) - 10.3(f) + 8.3(f) 

7 on 1 1 

^ “ 1^^^^ “ -\/31n(3)7r + —\/31n(3)27r + 0(e) . (B.2) 


The values for the remaining two-loop masters corresponding to M .^21 ^ M'S, MS and 

MS of Ref. [H] are trivial to construct as they correspond to products of one-loop masters, 

( 2 ) 

or the two-loop sunset graph in the case of ■ We note that the harmonic polylogarithms 
are based on the theory of cyclotomic polynomials [45] 
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